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I. INTRODUCTION 

A more general exponential screened Coulomb (MGESC) potential of the form: 

V{r) = - (^] [1 + (1 + br) exp(-26r)] , (1) 

where a is the strength coupling constant and b is the screening parameter, is known to 
describe adequately the effective interaction in many-body enviroment of a variety of fields 
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such as atomic, nuclear, solid-state, plasma physics and quantum field theory [1,2]. It is 
also used in describing the potential between an ionized impurity and an electron in a metal 
[3,4] or a semiconductor [5] and the electron-positron interaction in a positronium atom in 
a sohd [6]. 

The Schrodinger equation for such a potential does not admit exact solutions, various ap- 
proximate methods both numerical [7] and analytical [8] have been developed The MGESC 
Potential defined for an electron of the helium atom in the field of other electrons and nu- 
cleus has been investigated by Gerry and Laub [9]. Further, the large-N expansion was used 
to obtain the bound-state energy of the ground state and the first excited state and the 
corresponding wave functions analytically by Ref.[10]. 

In this paper, we calculate the binding energy eigenvalues of MGESC potential using a 
novel perturbative formalism [11] which has recently been used in solving the Schrodinger 
equation to obtain the bound-state energies as well as the wave functions for different types of 
potentials [11,12] in both bound and continuum regions. This novel treatment is based on the 
decomposition of the radial Schrodinger equation into two pieces having an exactly solvable 
part with an addiitional piece leading to either a closed analytical solution or approximate 
treatment depending on the nature of the perturbed potential. 

The contents of this paper is as follows. In Section II we breifly outline the method 
with all necessary formulae to perform the current calculations. In Section III we apply this 
method to the Schrodinger equation using a more general exponential screening Coulomb 
potential to obtain analytical expressions for the bound-state energy and eigen functions of 
different energy states. Finally, we end with some results and conclusions in Section IV. 

II. THE METHOD 

For a spherically symmetric potential, the corresponding Schrodinger equation, in the 
bound state domain, with the radial wave function reads 

^' <(^) - V(r) - E (2) 



with 



V{r) 



Voir) + 



K' 
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+ A\/(r), 



(3) 



2m r^ 

where AV(r) is a perturbing potential and ipn{r) = Xn{r)un{r) is the full radial wave func- 
tion, in which Xn{r) is the known normalized wave function of the unperturbed Schrodinger 
equation whereas M„(r) is a moderating wave function corresponding to the perturbing po- 
tential. Following the method given in Refs. [11,12], we may rewrite (2) as 
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The logarithmic derivatives of the unperturbed Xn(^) ^^'^ perturbed 'U„(r) wave functions 

are given by 
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where e„ is the eigenvalue for the exactly solvable potential of interest, and 
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-.AW:^ir) + 2Wnir)AW^ir) = A\/(r) - A£„, 



(7) 



in which Ae„ = E^p + E^^ + ■ ■ ■ is the correction term to the energy due to AVir) and 
En = £n + Asn- If Eq. (7), which is the most significant piece of the present formalism, can 
be solved analytically as in (6), then the whole problem, in Eq. (2) reduces to the following 
form 

h 



[Wnir) + AWnir) 



'2m, 



iWnir) + AWnir)y = Vir)-E„ 



(8) 



which is a well known treatment within the frame of supersymmetric quantum theory 
(SSQT) [13]. Thus, if the whole spectrum and corresponding eigenfunctions of the un- 
perturbed interaction potential are known, then one can easily calculate the required super- 
potential Wnir) for any state of interest leading to direct computation of related corrections 
to the unperturbed energy and wave function. 



For the perturbation technique, we can spht the given potential in Eq.(2) into two parts. 
The main part corresponds to a shape invariant potential, Eq. (6), for which the super- 
potential is known analytically and the remaining part is treated as a perturbation, Eq. 
(7). Therefore, the MGESC potential can be solved using this method. In this case, the 
zeroth-order term corresponds to the Coulomb potential while higher-order terms consitute 
the perturbation. However, the perturbation term in its present form cannot be solved ex- 
actly through Eq. (7). Thus, one should expand the functions related to the perturbation 
in terms of the perturbation parameter b, 

oo oo oo 

AVir;b) =J2bMr), AW-„(r;6) =J2b,W^\r), E^\b) =Y.b^E^\ (9) 

j=l i=l j=l 

where i denotes the perturbation order. Substitution of the above expansions into Eq. (7) 
and equating terms with the same power of b on both sides up to 0(6"^) yields 

2W^{r)Wi'\r) - ^^^f^"^ = V,{r) - E«, (10) 

v2m dr 



Wi'^^r) + 2W^{r)Wi'\r) - ^^^f^"^ = V,ir) - E^\ (11) 

V2m dr 

W^{r)W^\r) + Wi'\r)W^\r)\ - ^^^£^"^ = V,{r) - E^\ (12) 

Hence, unlike the other perturbation theories, Eq. (7) and its expansion, Eqs. (10)-(12), 
give a flexibility for the easy calculations of the perturbative corrections to energy and 
wave functions for the nth state of interest through an appropriately chosen perturbed 
superpotential. 

III. APPLICATION TO THE MGESC POTENTIAL 

Considering the recent interest in various power-law potentials in the literature, we work 
throughout the article within the frame of low screening parameter b. In this regard, the 
MGESC potential can be expanded in power series of the screening parameter b as [10,12] 



V{r) = - ( - ) [1 + (1 + br) exp(-26r)] = -- - - ^ V,ibr)\ (13) 

where the perturbation coefficients Vi are given by 

V, = -1, V2 = 0,Vs = 2/3, V, = -4/6, V, = 12/30, ■ ■ ■ . (14) 

Therefore, we apply this approximation method to the MGESC potential with the angular 
momentum barrier 
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where the ffist piece is the shape invariant zeroth-order which is an exactly solvable piece 
corresponding to the unperturbed Coulomb potential with Vo{r) = — 2a/r while AV{r) = 
ah — (2a6^/3)r^ + (4a6^/6)r^ — (12a6^/30)r^ + ■ ■ ■ is the perturbation term. The literature is 
rich with examples of particular solutions for such power-law potentials employed in different 
fields of physics, for recent applications see Refs. [14,15]. At this stage one may wonder why 
the series expansion is truncated at a lower order. This can be understood as follows. It is 
widely appreciated that convergence is not an important or even desirable property for series 
approximations in physical problems. Specifically, a slowly convergent approximation which 
requires many terms to achieve reasonable accuracy is much less valuable than the divergent 
series which gives accurate answers in a few terms. This is clearly the case for the MGESC 
problem [16]. However, it is worthwhile to note that the main contributions come from the 
ffist three terms. Thereby, the present calculations are performed upto the second-order 
involving only these additional potential terms, which suprisingly provide highly accurate 
results for small screening parameter h. 

A. Ground State Calculations (n = 0) 

In the light of Eq. (6), the zeroth-order calculations leading to exact solutions can be 
carried out readily by setting the ground-state superpotential and the unperturbed exact 
energy as 



Vr„=o (r) 
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(16) 



and from the literature, the corresponding normahzed Coulomb bound-state wave function 

[17] 
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X™ is an associate Laguarre polynomial function [18] 



(n— m)!(»n+fc)!r?i!' 

For the calculation of corrections to the zeroth-order energy and wave function, one needs 
to consider the expressions leading to the first- and second-order perturbation given by Eqs. 
(10)-(12). Multiplication of each term in these equations by Xnl^); ^^^ bearing in mind 
the superpotentials given in Eq. (5), one can obtain the straightforward expressions for the 
first-order correction to the energy and its superpotential: 
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and also for the second-order correction and its superpotential 
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for any state of interest. The above expressions calculate W'^\r) and W'^^\r) explicitly 
from the energy corrections E^^ and E^"* respectively, which are in turn used to calculate 
the moderating wave function Un{r). 

Thus, through the use of Eqs. (18) and (19), after some lengthy and tedious integrals, 
we find the zeeroth order energy shift and their moderating superpotentials as 
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Therefore, setting /9 = 6/a, the analytical expression for the ground s-state energy is explic- 
itly given, to order f3^, in atomic units {h = m, = 1) : 



Eja^ = _2 + /3 - i/3=^ + -/3^ - —^^ + 
0/ ^'^ 2*^ 8*^ 192'^ 



(22) 



and the full radial wavefunction is given by 

v2m 
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Hence, the explicit form of the full wave function in (23) for the ground state is 
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with P(r) = I]f=i Pi^* is a polynomial of fifth order having the following coefficients: 
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where a4 = 02 + ^w^ 
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B. Excited state calculations (n > 1) 

The calculations lead to a handy recursion relations in the case of ground states, however 
it becomes extremely cumbersome in the description of radial excitations when nodes of 
wavefunctions are taken into account, in particular during the higher order calculations. 
Although several attempts have been made to bypass this difficulty and improve calculations 
in dealing with excited states, (cf. e.g. [19], and the references therein) within the frame of 
SSQM. 

Using Eqs. (5) and (17), the superpotential Wn{r) which is related to the excited states 

can be readily calculated by means of Eqs. (18) and (19). Hence, the first-order corrections 

in the first excited state (n = 1) are 

^(1) ^ n'{i + 2f{i + 7){2i + 3) ^, 
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Consequently, substitution of the last equation into Eq. (19) allows us to write down 

48a'^m^ 
^^° {£ + 2f {£ + 3) {2£ + 3) {7f + 10^ + 211) ^^ 



ll52a*m^ 

Therefore, the analytical expressions for the first-excited s-state energy is explicitly given, 
to order (3^, in atomic units: 

E,/a' = -l+P-7P'+^-^P^-^P' + .... (28) 

The related radial wavefunction can be expressed in an analytical form by means of Eqs 
(18), (19) and (23), if required. The appromation used in this work would not affect con- 
siderably the sensitivity of the calculations. On the other hand, it is found analytically that 
our investigations put forward an interesting hierarchy between VF^^-* (r) terms of different 
quantum states in the first order after circumventing the nodal difficulties elegantly, 
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which, for the second excited state (n = 2) leads to the first-order correction 
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Hence, substituting W2 (r) into Eq.(19) gives the energy correction in the second-order as 
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Therefore, the analytical expressions for the second-excited s— state energy, to order f3^, in 
atomic units: 



F/n'- Kb ^^/?3+2295 347571 



(32) 



IV. RESULTS AND CONCLUSIONS 



Some numerical values of the perturbed energies of the Is, 2s and 3s states, in the atomic 
units, for various values of /3 in the range < /3 < 1.0 are presented in Table 1. Our results 
are consistent to order (3^ with earlier results obtained by applying the large A^-expansion 
method [10]. Therefore, our results are found in high agreement for small values of j3 with 
those given in [10]. On the other hand, for large screening parameter values, the accuracy of 
our results exceeds the ones given before in Ref.[10]. These results also appear to be in close 
agreement with the results obtained by solving Schrodinger equation numerically with the 
same potential via Numerov's method [10]. However, these results tend to deviate slightly 
as (3 approaches 1.0 with the results obtained from the numerical solution of Schrodinger 
equation. Moreover, we illustrate the improvement of energy with respect to orders of f3 in 
Table 2. 



In this work, we have shown that the bound-state energies of a MGESC potential for 
all eigenstates can be accurately determined within the framework of a novel treatment. 
Avoiding the disadvantages of the standard non-relativistic perturbation theories, the ob- 
tained expressions in the present work have much simple forms than the ones shown in the 
previous work [10]. 

Finally, the application of the present technique to MGESC potential for the first time is 
really of great interest providing leading to analytical expressions for both energy eigenvalues 
and wave functions which are likely of much interest in different field of physics. 
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TABLES 



TABLE I. Calculated binding energy eigenvalues for < /? < 1.0 up to order (3^. 
13 Eqq/o^'' Eoo/a"^ [10]^ Numerical [10] E^io/a^^ ^lo/a^ [10]^ Numerical [10] E20/C 



2a 



0.01 -1.9900005 

0.02 -1.9800039 -1.98000 (-1.98000) 

0.03 -1.9700130 

0.04 -1.9600304 -1.96003 (-1.96003) 

0.05 -1.9500586 -1.95006 

0.06 -1.9401161 -1.94010 (-1.94010) 

0.07 -1.9301565 

0.08 -1.9202305 -1.92023 (-1.92023) 

0.09 -1.9103236 

0.10 -1.9004377 -1.90044 (-1.90044) 

0.20 -1.8030143 -1.803 

0.30 -1.7086008 -1.70844 (-1.70958) 

0.40 -1.6169173 

0.50 -1.5269368 -1.523 (-1.537) 

0.60 -1.4374490 

0.70 -1.3477860 -1.321 (-1.384) 

0.80 -1.2587093 -1.2 

0.90 -1.1734581 

1.00 -1.0989583 -0.875 (-1.194) 



-0.4900067 
-0.4800516 -0.48005 
-0.4701668 
-0.460378 -0.46038 
-0.4507047 -0.45070 
-0.4411605 -0.44116 
-0.4317531 
-0.4224852 -0.42246 
-0.4133547 
-0.4043555 -0.4043 
-0.318752 -0.312 
-0.2431595 . 



(-0.48000) 

(0.46033) 

(0.44115) 

(-0.4221) 

(-0.4048) 

(-0.274) 



-0.2122538 
-0.2024526 

-0.183718 

-0.1662097 

-0.1495594 

-0.1334655 



^Calculations to order (3^. 
^Calculations to order /3^. 
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TABLE II. 


Improvement in 


energy with respect to orders of h/a. 




/3 


(^oo/«')o 


(i?oo/a2)i 


(.Eoo/a^)3 


{Eoo/a% 


{Eoo/a^h 


0.02 


-2.0 


-1.98 


-1.980004 


-1.9800039 


-1.9800039 


0.05 


-2.0 


-1.95 


-1.9500625 


-1.9500586 


-1.9500586 


0.08 


-2.0 


-1.92 


-1.920256 


-1.9202304 


-1.9202305 


0.20 


-2.0 


-1.80 


-1.804 


-1.803 


-1.8030143 


0.50 


-2.0 


-1.50 


-1.5625 


-1.5234375 


-1.5269368 


0.80 


-2.0 


-1.20 


-1.456 


-1.20 


-1.2587093 


(3 


{Eio/a^)o 


(£^io/a2)i 


{Eio/a% 


iEio/a^)4 


{Eio/a^)6 


0.02 


-0.5 


-0.48 


-0.480056 


-0.4800516 


-0.4800516 


0.04 


-0.5 


-0.46 


-0.460448 


-0.4603776 


-0.460378 


0.05 


-0.5 


-0.45 


-0.450875 


-0.4507031 


-0.4507047 


0.06 


-0.5 


-0.44 


-0.441512 


-0.4411556 


-0.4411605 


0.08 


-0.5 


-0.42 


-0.423584 


-0.4224576 


-0.4224852 


0.10 


-0.5 


-0.40 


-0.407 


-0.40425 


-0.4043555 


0.20 


-0.5 


-0.30 


-0.356 


-0.312 


-0.318752 


0.30 


-0.5 


-0.20 


-0.389 


-0.16625 


-0.2431595 


/? 


(^2o/a^)o 


{E2o/a\ 


{E2o/a^)2 


iE2o/a^)4 


iE2o/a^)6 


0.02 


-0.2222222 


-0.2022222 


-0.2024982 


-0.2024523 


-0.2024526 


0.04 


-0.2222222 


-0.1822222 


-0.1844302 


-0.1836958 


-0.1837180 


0.06 


-0.2222222 


-0.1622222 


-0.1696742 


-0.1659563 


-0.1662097 


0.08 


-0.2222222 


-0.1422222 


-0.1598862 


-0.1481358 


-0.1495594 


0.10 


-0.2222222 


-0.1222222 


-0.1567222 


-0.1280347 


-0.1334655 
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